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We present lattice results for the gluon propagator for SU(2) and SU(3) in the 
CIh' Laplacian gauge which avoids lattice Gribov copies. In SU(3) we compare with 

the most recent lattice calculation in Landau gauge and with various approximate 
CaJ ■ solutions of the Dyson Schwinger equations (DSE). 
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Introduction 



f\J ■ We first summarize the results obtained within the Landau gaug«l: By 

solving approximately the DSE, Mandelstam found an infrared enhanced gluon 

f— n | propagator of the form D(q 2 ) ~ -j-. Avoiding gauge copies, Gribov obtained 

^2 \ D(q 2 ) ~ i+ m i ■ Using the "pinch technique", Cornwalla obtained a solution 

which fulfills the Ward identities, allows a dynamical mass generation, and also 
predicts a finite value for D(0) = D(q 2 — 0) consistent with our data. 



Early results for the gluon propagator obtained directly from Lattice QCD 
on small lattices u were interpreted in terms of a massive scalar propagator. 
Results on larger lattices were accounted for by assuming a positive anoma- 
Q_i! lous dimensionH: D(q 2 ) ~ 2(1+ ^ 2 . A recent, detailed study of the gluon 

propagator uses very large lattices Q. Since we want to compare our results 
with these, we follow closely their analysis and refer to Refs. Bu for details. 
In the Laplacian gauge, the longitudinal part of the gluon propagator does 
not vanish; the transverse scalar function D(q 2 ) can be extracted from T>°^(q) 

as D(q 2 ) = Ue^|E«^(«)} " l^fi > where F ^ 2 ) is determined by 
projecting the longitudinal part of T>"^(q) using the symmetric tensor q^q v . 
Gauge Fixing Procedure 

Previous lattice studies all fixed to Landau gauge by using a local iterative 
maximization algorithm, which converges to any one of many local maxima 
(lattice Gribov copy), but fails to determine the global one. To overcome 
this problem, we use a different gauge condition, the Laplacian gaugea, which 
is Lorentz-symmetric and gives a smooth gauge field like the Landau gauge, 
but which specifies the gauge unambiguously. We consider the maximization 
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of Q = Re^ x Tr \g{x)U ll {x)g J{ (x + fi) — g(x)g(xY] . If one relaxes the re- 
quirement that g S SU(N), maximizing Q is equivalent to minimizing the 
quadratic form J2 xy fx^xyfy, with A(U) the covariant Laplacian. Using the 
(N — 1) lowest-lying eigenvectors fi(x) of A(U), one can fix the gauge uniquely 
by requiring Vx,/|(x) eR,//(jc) = 0, j = (i + 1), .., N%. 
Results 

In Fig.l we show the transverse gluon propagator for SU(2) Yang-Mills 
theory in two different volumes; too = ^-D(O) -1 for the 16 4 lattice. Changing 
the volume has little effect, in particular on D(0). We observe similarly small 
volume effects in SU(3). This is strikingly different from Landau gauge, where 
Zwarrziger has argued that -D(O) should vanish in the infinite lattice volume 
limit a. This prediction is indeed consistent with recent lattice results in SU(2) 
at finite temperatureH. In contrast, in the Laplacian gauge, we find that D(0) 
is finite and independent of the volume V for V larger than about l/2fm ~ 
D{0) 2 . We find D(0) = 58(2) in lattice units at (3 = 6.0, i.e. D(O)" 1 / 2 = 248(5) 
MeV (using a -1 = 1.885 GeV), corresponding to a length scale of about 0.8 fm. 

In Fig. 2 we compare results for the gluon propagator in SU(3) quenched 
QCD in Laplacian and Landau gauges, (too = a/D(0) _1 in the Laplacian 
gauge). Scaling is checked on the 16 3 x 32 lattice for /3 = 5.8 and 6.0. Making 
a cylindrical cut in the moment a Q to minimize lattice artifacts, we find that 
scaling is very well satisfied for the Laplacian gauge, with both sets of data 
falling on a universal curveQ. 

We fit to our data the same models as considered by Leinweber et al. U 
in Landau gauge. Since we have observed scaling, we use our results at the 
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Figure 1: The SU(2) gluon propagator in 

two different volumes. 
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Figure 2: The SU(3) gluon propagator in 
Laplacian and Landau gauges. 
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D(0) 


X 2 /d.o.f 


Gribov 


2.63(2) 


0.203(7) 
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Stingl 


2.63(2) 


0.203(13) 


0.002 (1.100) 
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Marenzoni 


2.47(3) 


0.199(6) 


0.237(5) 
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Cornwall 


7.08(9) 


0.281(4) 


0.265(8) 




59 


2.5 


Model A 


1.96(1) 


0.654(17) 


2.181(67) 


8.91(41) 


43 


1.2 



Table 1: best fit of parameter values to our /3 = 6.0 data on the 16 3 x 32 lattice. 

finer lattice spacing (f3 = 6.0) for the fits. Table 1 and Fig. 3 summarize the 
results of the fits to the various models. We find that Gribov-type models are 
excluded, whereas Cornwall's model is clearly favored among all analytically 
motivated models. Model "A" Q, which gives a better fit, is phenomenological, 
contains one more parameter, and misses D(0) by 25%. One can then use the 
fit to Cornwall's model to analytically continue to negative q 2 and determine 
the gluon pole mass. This is carried out in Ref. EL 
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In conclusion, we see significant 
modifications from Landau gauge in 
the infrared. In particular, we find 
that -D(O) obeys scaling, is finite, 
and volume independent for large 
enough volumes. We find support 
for Cornwall's model which fits the 
momentum dependence of the prop- 
agator rather well, whereas models 
with infrared enhancement of the 
type l/(g 2 ) 2 or Gribov-type sup- 
pression are excluded. 
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